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ABSTRACT 

We examine the incidence of cold fronts in a large sample of galaxy clusters extracted from a 
(512/i _1 Mpc) hydrodynamic/N-body cosmological simulation with adiabatic gas physics computed 
with the Enzo adaptive mesh refinement code. This simulation contains a sample of roughly 4000 
galaxy clusters with M> 1O 14 M at z=0. For each simulated galaxy cluster, we have created mock 
0.3-8.0 keV X-ray observations and spectroscopic-like temperature maps. We have searched these 
maps with a new automated algorithm to identify the presence of cold fronts in projection. Using 
a threshold of a minimum of 10 cold front pixels in our images, corresponding to a total comoving 
length L c f > 156h~ 1 kpc, we find that roughly 10-12% of all projections in a mass-limited sample 
would be classified as cold front clusters. Interestingly, the fraction of clusters with extended cold 
front features in our synthetic maps of a mass-limited sample trends only weakly with redshift out 
to z=1.0. However, when using different selection functions, including a simulated flux limit, the 
trending with redshift changes significantly. The likelihood of finding cold fronts in the simulated 
clusters in our sample is a strong function of cluster mass. In clusters with M > 7.5 x 10 14 M Q the 
cold front fraction is 40-50%. We also show that the presence of cold fronts is strongly correlated 
with disturbed morphology as measured by quantitative structure measures. Finally, we find that the 
incidence of cold fronts in the simulated cluster images is strongly dependent on baryonic physics. 
Subject headings: Galaxies: clusters: intracluster medium-large-scale structure of Universe-X-rays: 
galaxies: clusters-Methods: numerical 



1. INTRODUCTION 

With the advent of high resolution X-ray imaging using 
space-based observatories, in particular Chandra, unex- 
pected discontinuous X-ray features in galaxy clusters 
were discovered. The expected appearance of shocks 
was supplemented by the appearance of cold fronts, 
features that do not exhibit a pressure jump, in con- 
trast to shocks, but have jumps in surface brightness 
and temperature in opposing directions (for review see 
iMarkevitch fc Vikhlinin! [2007). Temperature jump ra- 
tios (from one side of the feature to the other) in typ- 
ical cold fronts range from factors of 50% to factors of 
a few, with corresponding (but reversed) inferred den- 
sity jumps. Because of the pressure continuity across 
these features, they are typically described as contact 
discontinuities. Early classic cold fro nt examples in- 
clude Abell 3667 (VikhlininelaLl Hooi3) and Abell 2142 
([Markevitch et al.ll2000l ). The properties of cold fronts 
have turned out to provide excellent constraints on the 
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physics of the intracluster medium (ICM). For exam- 
ple, the extremely sharp edge (typically narrower than 
the Coulomb mean free path) between the cold and hot 
sides of the features places limits on the effectiveness of 
conduction acro ss the edge (e.g.. lEttori fc Fabianl f2000: 
lAsai et ah! 120041 ) . The apparent lack of fluid instabil- 
ity growth at the interface also suggests some physical 
process that stab ilizes the edge, perhaps the draping of 
magnetic fields (I Vikhlinin et al.l l2001bl : iLvutikovl 120061 : 
lAsai et al.ll2007l : lTakizawall2008T ). 

Following the early discoveries of cold fronts, many 
other clusters have been identified as hosting cold 



fronts (e.g.jMazzotta et a l. 2001; Markevit ch et all 2002: 
Dupkefc WMt3 120031 : [Hallman & Markevitchl l2QQ4t 

~ JboiflD- 



Johnson et al.l I2Q1QI ) . Statistics of cold fronts in clus- 
ters hav e been calculated u sing data from the Chandra 
archive (|Owers et al.l I2009D and from flux- l imite d sam - 
ples using XM M-Newton (iGhizzardi et all 120061 . l2Q10h 
and Chandra ([Markevitch et al.l l2QQ3h . The main re- 
sults are that in flux-limited samples, anywhere between 
40% and 87% of clusters ar e found to host cold fronts 
([Ghizzardi et al.l 120061 l2QlQh depending on the redshift 
range explored. Also, a lar ge fraction of cool core c lusters 
(-67%) host cold fronts (IMarkevitch et all 120031 ). Re- 
cent work also suggests that the appearance of cold fronts 
is strongly correlated with cluster mergers ([Owers et al.l 
120091 : iGhizzardi et al.ll20Toh . 

Numerical hydrod ynamic simulations of i d ealized clus- 
ter m ergers (e. g .. iRoettiger et al.l Il999l: iHeinz et Id] 



20031: iTakizawal l2QQ5t j Ascasibar & Markevitch 2006; 
Mastropietro fc Burke rt 2008; Soring el fc Farrarl 120071 
ZuHone fc M arkevitch 2009) as well as cosmological sim- 
ulations (e.g., Nagai fc Kravtsovll2003l : lBialek et al .11200 2l 



iMathis et aD2005l iTittlev fc Henrikseril2005l ) appear to 
reproduce the properties of all types of observed cold 
fronts naturally in the process of mergers. Recent stud- 
ies have begun to separately classify two types of cold 
fronts, the first being "merger" cold fronts, where the 
system is obviously disturbed in X-ray morphology, and 
clearly undergoing a c urrent major merger, an d so-called 
"sloshing" cold fronts (M arkevitch et al.ll200ll ) in appar- 
ently otherwise relaxed cool-core clusters. This distinc- 
tion is somewhat artificial, as simulations show that both 
cold front types are associated with mergers, or in the 
case of sloshing, at least a close pass of a subcluster. 

Simulation studies indicate that cold fronts arise from 
one of a few physical scenarios. First, and most obvious 
are the so-called merger- type cold fronts. In these sce- 
narios, two objects merge, at least one a massive cluster. 
In the early stages, while the merging subclusters are 
approaching, their relative velocity is supersonic (in the 
cluster gas) and drives shocks. In the simplest scenario, 
the ram pressure associated with the relative velocity 
of the subcluster in the gas of the main cluster pushes 
the subcluster gas out of the dark matter potential well. 
Once the gas separates from the dark matter potential, 
it expands adiabatically and cools, resulting in the cold 
front feature. In this scenario, cold fronts should always 
be associated with shocks, which may also be visible in 
the ICM adjacent to the cold front. Depending on when 
in the process of merging the cold front is observed, the 
relative positions of cold gas, dark matter, and shocks 
can vary. 

The sloshing type cold fronts are observed typically 
in cool core clusters, and appear at small cluster radius 
(R< lOOkpc). These clusters otherwise appear relaxed 
dynamically, leading to the separate classification from 
merger type cold fronts. There is no obvious evidence of 
merging in these clusters in the X-ray images, but simu- 
lations that reproduce cold fronts of this type show that 
they are also caused by dynamical interactions. In simu- 
lations, a near pass of a subcluster pulls the cluster mass 
away from the original center position, and as the sub- 
cluster moves away, the central baryonic material rocks 
back and forth in the dark matter potential. This creates 
a characteristic spiral shaped cold front near the center 
of the cool core cluster. 

Both the simulation work referenced here and recent 
theoretical study (e.g., jLvutikovl 120061 : iBirnboim et al.l 
■Ll l2009h 



2009; iKeshet et al.ll2009f ) have made significant progress 
toward understanding the process by which cold fronts 
form in galaxy clusters and the detailed physics that de- 
termines their observed properties. What has not been 
attempted to date is an estimate of the predicted inci- 
dence of cold fronts in clusters from fully cosmological nu- 
merical hydro/N-body simulations with large (N > 1000) 
samples. The appearance of merger- type cold fronts is 
governed by both the merger rate of cluster scale systems, 
as well as the gas dynamics and more subtle effects like 
projection angle. Therefore, in order to determine how 
many clusters we expect to have cold fronts at a given 
epoch, or mass limit, or depth of observation, a statistical 
treatment is required. Though idealized merger scenar- 
ios, or small samples simulated at high resolution are 
critical to understanding the range of cold front forma- 
tion, they neglect to examine the full ranger of mergers in 
large scale structure. Fully cosmological simulations, like 



the ones performed for this study, create realistic merger 
histories, including multiple mergers. From these simu- 
lations we can gain valuable insight about the process of 
cold front formation and evolution over a large statistical 
sample of galaxy clusters. 

The outline of this paper is as follows. In Section [2j we 
describe the setup and analysis of the Enzo simulations. 
In Section [3] we discuss the method of cold front identi- 
fication. In Section |4] we describe our major results. In 
Section [5] we describe briefly the effects of numerical res- 
olution and baryonic physics. Finally, in Section [6j we 
discuss and summarize our results. 

2. SIMULATION SETUP AND ANALYSIS 

The main statistical results come from the so- 
called Sa nta Fe Light Cone ( S FLC) simulation, de- 
scrib ed in lHallman et afl (120071) . iSkillman etHI (|2008h 
and lHallman et al.l |2009l ). This calculation is per- 
formed with 'Enzo , a publicly available adaptive 
mesh refinement (AMR) cos mology code developed b y 



mesn rennement (AMU) cosmolo gy code developed Dy 
Greg Bryan and colleagues ([B ryan & Norman 1997aTlbl; 
iNorman fc Brvar] [l999l : lO'Shea et al.ll2004 l2005h . The 
specifics of the Enzo code are described in detail in these 
papers (and references therein). 

The E nzo code couples a n N-b o dy particle- mesh (PM) 
solve r ([Efstathiou et al.l 1 19851 : iHocknev fc Eastwood! 
1988) used to follow the evolution of a collisionless dark 
matter component with an Eulerian AMR met hod for 
ideal gas dynamics by iBerger fc Colellal (|1989[ ). which 
allows high dynamic range in gravitational physics and 
hydrodynamics in an expanding universe. 

This simulation is set up as follows. We initialize our 
calculation at z = 99 assuming a cosmological model 
with Q m = 0.3, Q b 0.04, ^ A 0.7, h 0.7 
(in units of 100 km/s /Mpc), as = 0.9, and using an 
lEisenstein fc Hul (|l999l ) power spectrum with a spectral 
index of n s — 1. The simulation is of a volume of the uni- 
verse 512 h _1 Mpc (comoving) on a side with a 512 3 root 
grid. The dark matter particle mass is 7.3 x 10 10 h _1 M . 
The simulation was then evolved to z = with a max- 
imum of 7 levels of adaptive mesh refinement (a maxi- 
mum spatial resolution of 7.8 h _1 comoving kpc), refin- 
ing on dark matter and baryon over densities of 8.0. The 
equations of hydrodynamics were solved with the Piece- 
wise Parabolic Method (PPM) using the dual energy for- 
malism. This simulation is performed using adiabatic 
physics only. It has become clear from observations that 
additional non-gravitational physics is important in the 
evolution of the ICM. However, this run models a large 
volume with high peak resolution, which is only now be- 
coming possible to calculate while including additional 
baryonic physics. It has the advantage of generating a 
cluster sample with thousands of massive clusters, allow- 
ing us to explore a wide range of cluster interactions that 
result in cold fronts. 

Analysis was performed on 11 data outputs between 
z = 1 and z = 0.0 (spaced with Sz = 0. 1) in an identical 
way. The HOP halo-finding algorithm (Ei senstein fc Hutl 
1998) was applied to the dark matter particle distri- 
bution to produce a dark matter halo catalog (as in 
lHallman et al.ll2007f ). For each output, we calculate bulk 
properties (e.g., total mass, X-ray luminosity), radial 

8 http: / /lca.ucsd.edu /projects / enzo 
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Fig. 1. — Histogram of number of clusters with M > 1O 14 M0 in 
the comoving volume as a function of redshift. Shows the sample 
size used for the cold fronts analysis. 

profiles, and projected images (along 3 orthogonal axes) 
of each simulated cluster in a variety of physical and ob- 
servable properties. Spherically-averaged, mass- weighted 
radial profiles of various baryonic and dark matter quan- 
tities including density, temperature, and pressure were 
generated for every halo in the catalog with an estimated 
halo mass greater than 4 x 10 13 M . These radial pro- 
files were used to calculate M200 and i?200- ^200 refers 
to the total mass inside a radius of i?200, the radius at 
which the overdensity average inside the sphere centered 
on the cluster is 200 times the critical density. We then 
generate synthetic X-ray images for clusters in the sam- 
ple with M200 > 1O 14 M . We then run our automated 
cold front finder on three orthogonal projections of each 
cluster. A histogram of the size of the mass-limited clus- 
ter sample as a function of redshift is shown in Figure [TJ 

2.1. Testing the Effect of Baryonic Physics 

For the purpose of understanding the impact of addi- 
tional baryonic physics on the cold front incidence, we 
have analyzed two additional simulations. Each of the 
two simulations have identical simulation volumes, ini- 
tial conditions, and cosmology. The simulated volume 
is a 128/z -1 comoving Mpc cube, with 256 3 root grid 
zones and 256 3 dark matter particles. Each run is al- 
lowed to refine to a maximum of 5 AMR levels, and the 
refinement criteria are overdensity of 8.0 in either the 
dark matter or the gas (with respect to the mean on the 
parent level). This gives the run a peak spatial resolu- 
tion of 15.6/i _1 kpc. The dark matter particle mass is 
9.0 x 10 9 h _1 M and the mean baryon mass resolution 
is 1.4 x 10 9 h _1 Mq, a factor of 8 improvement over the 
SFLC run used for our main analysis. For these runs, we 
use a cosmological model with Q m = 0.27, ^5 = 0.044, 
Q A = 0.73, h = 0.7 (in units of 100 km/s/Mpc), a 8 = 0.9, 
and using a power spectrum with a spectral index of 
n s = 0.97. These values are changes from our previ- 
ous run, to match WMAP5 cosmological parameters. 
Therefore, there are multiple changes to these simula- 
tions from the SFLC, including the mass and spatial res- 



olution, baryonic physics, and cosmological parameters. 
The most relevant changes to the cosmology are a slight 
reduction in the value of ^ m and a slight increase in 
the value of However, in this work, we endeavor to 
merely compare our two runs where everything is identi- 
cal except baryonic physics, to examine the impact. 

The two additional simulations presented here are 
identical in terms of the initial conditions and simulation 
properties, except that one has only adiabatic and gravi- 
tational physics, and the other includes a prescription for 
radiative cooling of the gas using non-eq uilibrium cool- 
ing a nd chemistry for H and He, Cloudy (jFerland et al.l 
1998) for the metal cooling, and star formation includ- 
ing ther mal and metal fe e dback from supernovae as de- 
tailed in ICen fc Ostrikerl (|1999[ ). The cooling plus star 
formation and feeback run includes a spatially uniform 
but time- varying ultravi olet (UV) radiation background 
(jHaardt fc Madaul fl996h . We turn on the UV back- 
ground at a redshift z = 7. For the metal line cooling, 
we interpo- late from a grid of data made with Cloudy. 
For redshifts z > 7, we use a similar grid of heating 
and cooling data, without the influence of the UV ra- 
diation background, assuming collisional ion ization only 
(for more details, see ISmith et al.ll2008Ll201Q[ ). Addition- 
ally, we have highly time-resolved data outputs for these 
calculations, allowing us to examine the time history and 
origin of cold fronts in the simulated clusters in detail. 
In this paper, we examine only the general incidence of 
cold fronts in these simulations, and leave a study of the 
detailed properties of the cold fronts, as well as a full 
convergence study, to future work. 

For our simulation with purely adiabatic physics (CC- 
Adia), and the identical run with additional baryonic 
physics (CC-OTH), the data outputs and analysis are 
identical. We perform the same analyses as done for the 
SFLC simulation. 

3. COLD FRONT IDENTIFICATION AND ANALYSIS 

The identification of X-ray cold fronts in the syntheti- 
cally observed numerical clusters is done via a relatively 
simple algorithm, which is designed to roughly match 
the observational (by eye) identification, but in an auto- 
mated way. 

For this study, we focus on the identification of X- 
ray cold fronts, characterized as features which are both 
colder and brighter than the surrounding intracluster 
medium. This identification is made from projected im- 
ages of the cl usters of the spectroscopic-like tempera - 
ture (T 3 i) (see iMazzotta et all 12003 : iRasia et aTll20Q5D . 
and the 0.3-8.0 keV X-ray surface brightness calcu- 
lated using the Cloudy code (Fer land et afl [l998). The 
spectroscopic-like temperature has been determined to 
be an accura te proxy for the m easured X-ray spectral 
temperature (Ra sia et al.l l2005h . We are not using the 
vo lume integra l of th e spectroscopic-like temperature as 
in IRasia et aD (|2005f ). but are using the same weight- 
ing for the line of sight integral in order to make mock 
temperature maps. This weighting has been shown to 
reproduce the fitted spectral temperature better than 
standard emission weighting or mass weighting. The cal- 
culation of T s i for our projected maps is 

_ Jn 2 T a /T°- 5 dl 
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where a=0.75. 

For each image, we step in each row or column of the 
pixel map, and look for jumps in the surface brightness 
and temperature which are above a threshold value, and 
are opposite in sign. The jump thresholds are 

S Xi i/S Xt2 > 2.0, (2) 
T sh2 /T shl > 1.4, (3) 

where the numerical subscripts indicate the 2 opposing 
sides of the edge in temperature and surface brightness. 
Note that the jump in surface brightness is in the oppo- 
site direction of the jump in temperature. Observed cold 
fronts have temperature jumps from factors of 50% to 
factors of a few, so this threshold should be sufficient to 
recover the observed distribution. The threshold jump 
strengths are designed to be roughly continuous in pres- 
sure, as the fixed band X-ray surface brightness is only 
weakly dependent on temperature and depends on the 
electron density squared. Therefore a surface brightness 
jump of 2.0 is roughly equivalent to a projected density 
jump of 1.4. We also limit the identification of cold fronts 
to the area within a radius of r$oo projected on the sky, 
where rsoo is the radius within which the mean overden- 
sity of the cluster is 500 times the value of p c , the critical 
density of the universe. This radius is chosen since it is 
representative of the part of massive galaxy clusters that 
can typically be observed with standard exposures using 
current X-ray telescopes. 

We expect the number of cold fronts and the size of 
the jumps to be a lower limit to the true number in 
the simulation. This results from a number of factors. 
First, the effects of projection on the detection of cold 
fronts must be considered. Like shocks in galaxy clusters, 
the detection is strongly orientation dependent. Edge-on 
orientations lead to much more detectable shocks and 
cold fronts, face-on oriented features are basically un- 
detectable via these methods. Second, because of both 
orientation and potential spreading of the discontinuity 
over several numerical grid cells, the jump may span 2-3 
pixels. Our method in that case may miss a more grad- 
ual jump, or detect it with a smaller ratio of temperature 
and surface brightness across the feature. The first effect 
(projection) is not considered directly in this study, since 
we only wish to compare the observed incidence of cold 
fronts to the simulated incidence, and real cluster obser- 
vations suffer the same effect. We partially compensate 
for this effect by using multiple orthogonal projections 
of each cluster, but this doesn't impact the fractional in- 
cidence, as we count each projection as an independent 
image. The second set of effects is more challenging to 
correct. In this study we have experimented with check- 
ing for features matching our criteria across both 2 and 
3 pixels. The main result of that experiment is that the 
overall statistics are not strongly dependent on the num- 
ber of pixels across which the jump is considered. This 
results from the fact that jumps that are identified in a 
single pixel jump are also typically identified in 2-3 pixels 
as well. For this same reason it is not effective to include 
all jumps of 1, 2, and 3 pixels for instance in our analysis 
due to redundant detections of the same features. For 
the purposes of this study, we consider the incidence of 
cold fronts detected via this method to be lower limits, 
and the jumps may also be underestimated (individu- 
ally). Based on our tests, we expect the overall statistics 



to be largely unaffected, though it will be explored in 
more detail in our later work. 

Because cold fronts are extended features in the X-ray, 
we anticipate that their incidence will not be strongly de- 
pendent on grid resolution. The incidence of cold fronts 
may however be somewhat dependent on the baryonic 
physics of the simulation, as well as the mass resolution. 
Mass resolution of the simulation results in a higher num- 
ber of resolved subclusters, increasing the rate of merg- 
ers. Although previous work has shown that cold fronts 
certainly appear in adiabatic simulations, purely adia- 
batic processes may not be the only mechanisms for cold 
front creation. The so-called "sloshing" cold fronts may 
result from gas which has radiatively cooled in the cen- 
ters of clusters oscillating in the dark matter potential in 
response to mergers. These sloshing features appear in 
obser ved clusters with strong central entropy grad ients 
(e.g., iMarkevitch et al.l [20031 : iGhizzardi et al.ll201Q[ ). In 
this work we make no distinction between sloshing and 
merger- type cold fronts, we merely identify all features 
that match the cold front criteria specified above. In our 
adiabatic simulations, because we do not get strong core 
entropy gradients (as are seen in the presence of strong 
radiative cooling), we expect our cold fronts to be of the 
merger type in this simulation. 

4. RESULTS 

4.1. Incidence of Cold Fronts in Simulated Clusters 

The first result of the cold front identification is that 
in the adiabatic simulation, for the sample of all clusters 
with M > 10 14 M o out to z = 1, 12% of images have 10 or 
more cold front pixels, 5% have 20 or more, and roughly 
1% have 50 or more. At z=0 in our adiabatic simulation, 
we have 116 images with more than 50 cold front pixels. 
In these images, there are clear, unambiguous cold fronts 
spanning large areas of the cluster projection (see Figure 
[2|). Figure [3] shows the result of counting clusters above 
some minimum pixel count identified as cold fronts. This 
is equivalent to measuring the comoving total length of 
cold front features in each image, as each pixel has a con- 
stant comoving size in all of these images. In the figure, 
we plot the trend as a function of redshift for clusters 
with at least 10 pixels (L c f > 156/z -1 kpc total) identi- 
fied as a cold front, 20 pixels (L c f > 312h~ 1 kpc) and 50 
pixels (L c f > 7S0h~ 1 kpc). These numbers are arbitrary, 
and designed purely to show the general redshift trend. 
It is clear from this plot that there is only a weak trend 
in the incidence of cold fronts as a function of redshift 
out to z =1. 

The choice of a threshold in cold front pixels or comov- 
ing length of features is fairly arbitrary, but identification 
of one pixel with cold front criteria is not particularly in- 
teresting. First, it is quite possible to have spurious cold 
pixels due to line-of-sight superpositions. Additionally, 
cold fronts are extended features, and so in these images 
at least several pixels need to be associated in order for 
a true identification to take place. 

Another interesting statistic is the incidence of cold 
fronts in simulated clusters as a function of mass. For 
this case, we take clusters with at least 10 identified cold 
front pixels, which as described earlier is equivalent to 
the combined length of cold fronts L c f > 156h~ 1 kpc. 
We combine the result for all clusters from 1.0 > z > 0.0 
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Fig. 2. — Representative cold fronts from the numerically 
simulated clusters at z=0, with mass (from top to bottom) of 
3.1xlO 15 M , 1.4xlO 15 M , and 1.2xlO 15 M . Image quantity is 
T s i, contours are 0.3-8.0 keV X-ray surface brightness. Field of 
view is roughly 3h _1 Mpc. Image scale goes from T ~3keV (blue) 
to T ~10-12keV (yellow- white). Arrows indicate locations identi- 
fied by the automated finder as cold fronts. 

in Figure 2J What we see here is that the likelihood 
of finding extended cold fronts in clusters increases with 
mass. There are several possible interpretations for this 
result. The first, which is that this is a real physical 
trend, is that more massive clusters form in more dense 
local environments, leading to more activity from merg- 
ers and accretion, increasing the number of cold fronts 
generated. A second possibility is that because smaller 
clusters are smaller in projection, and we have a fixed co- 
moving pixel scale, they are less likely to have extended 
cold front features purely because of resolution effects. 
This effect is compounded by the larger cross section for 
mergers for the larger clusters, generating a higher likeli- 
hood for a merger of a given size in that projected volume 
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Fig. 3. — Fraction of clusters displaying a cold front as a function 
of redshift. Solid line is clusters with more than 10 cold front pixels 
(total L c f > 156/i -1 kpc comoving), dotted line is for 20 cold front 
pixels (L c f > 312h~ 1 kpc), the dashed line represents 50 pixels 
(L c f > 780h~ 1 kpc) identified as cold fronts in the image. Gray 
regions are la Poisson error bars. Note the relatively flat redshift 
trend. 

at any time. It is legitimate to ask how the sample can 
have only a weak trend in cold front incidence with red- 
shift, yet have a strong mass trend. Given the hard lower 
mass limit at each redshift, the distribution of clusters 
as a function of mass in each redshift bin is different. In 
other words, the clusters selected at each redshift with 
an identical mass cutoff do not represent the same part 
of the halo mass function at each redshift. This effect 
may contribute to the result shown here. We expect that 
if cold fronts are related to merger activity (as is indi- 
cated by simulation studies) that their incidence should 
trend with redshift, as large-scale structure grows, and 
the merger rates change as a function of time. 

In our above calculation of the fraction of cold front 
clusters as a function of cluster mass, we include clus- 
ters from all redshifts. Does this trend with mass vary 
as a function of redshift? Figure [5] shows the fraction of 
cold front clusters as a function of mass for 4 redshift 
bins in the simulation. The result is that the trend with 
mass is nearly identical for all redshifts, albeit with slight 
shifting in the upper mass bins due to small numbers. Ef- 
fectively, clusters at all redshifts show an identical trend 
of increasing fraction of clusters with cold fronts with 
cluster mass. 

4.2. Selection Effects 

Using a mass or flux-selected sample of galaxy clusters 
to study trends as a function of redshift is an inherently 
flawed method, as has been pointed out by many investi- 
gators. The most well known bias in flux-limited samples 
is the Malmquist bias, where your survey progressively 
selects brighter objects at greater distances. For galaxy 
clusters this has the effect of choosing the higher mass 
objects at higher redshift. Our sample has a bias working 
in the opposite direction, which is that we are selecting a 
mass- limited sample from all redshifts. Effectively we are 
sampling objects in different parts of the halo mass func- 
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Fig. 4. — Trend of cold front incidence as a function of mass 
for clusters from 1.0 > z > 0.0. Solid line is for clusters counted 
as cold front clusters if they have at least 10 identified cold front 
pixels, corresponding to a total cold front length L c f > 156h~ 1 kpc 
comoving, within a projected radius of rsoo- These pixels are not 
necessarily contiguous. Also plotted for comparison are the trends 
with a limit of half as many pixels (L c f > 78h~ 1 kpc, dotted line) 
and twice as many (L c f > 312/z _1 fcpc, dashed line). Note the steep 
trend in mass, and the flattening at high mass (M > 7.5x 1O 14 M0). 
Gray regions are la Poisson error bars. 
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Fig. 5. — Fraction of clusters with L c f > 156h~ 1 kpc as a function 
of cluster mass for 4 different redshift bins, z=0, z=0. 3, z=0. 5 and 
z=0.8. Blue dashed lines are for la Poisson errors for the z=0 
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end. 



Poisson errors overlap for all redshifts at the high mass 



tion at each redshift, not objects that at high redshift 
would be the precursors to the objects at low redshift in 
our sample. The result of such a selection is that true 
evolutionary effects as a function of redshift can easily 
be masked . One solution is to use a method similar to 
lHart et all ([2009) , selecting objects along an evolution- 
ary "road" , such that objects selected at low redshift are, 
in a statistical sense, the children of the objects at higher 



To understand our result of the incidence of cold 
fronts having only weak trending with redshift, we sort 
our cluster sample in two additional ways. The first 
method we choose here is to use a Press-Schechter 
([Press fc Schechterlll974l ) mass function analysis to se- 
lect a lower mass cutoff at each redshift corresponding 
to halos of the same rarity in large scale structure. This 
threshold choice effectively follows the same cluster pop- 
ulation as a function of redshift, allowing us to see evo- 
lutionary trends in the incidence of cold fronts. 

Our selection is done by calculating the value of (Jm, 
the RMS density fluctuation, as a function of cluster 
mass and redshift. We then choose a lower mass limit for 
our sample such that gm is a constant with redshift. This 
choice is motivated by the analytic expressions for the 
cluster mass function, both derived from linear theory 
and fit empirically to numerical simulations. The deriva- 
tions below are taken from the cited references as well as 
lHallman et al.l (|2QQ7f ). The comoving numbe r density of 
clusters as expressed bv I Jenkins et al.l (|200l[ ) is 



dn 
dm 



(M, z) = -0.315^ — ^ 
v } M cfm dm 

exp 



[_| .61-^P(^m)| 3 - 8 ] (4) 

where gm is the RMS density fluctuation, computed on 
mass scale M from th e z = linear power spectrum 
([Eisenstein fc HuHl999[ ). po is the mean matter density 
of the universe, defined as po = ft m pc (with p c being the 
cosmological critical density, defined as p c = 3Hq / 8ttG) , 
and D(z) is the linear growth function, given by this 
fitting function: 



D(z) 



tor. 



_ 1 5£l m (z) 
~ 1 + z 2 

{zf 7 - n A (z) + [i 



n m (z) o.a(z)- 



70 



(5) 



([Carroll et al.l Il992[ ). with Q m (z) and ^a(z) defined 
in the typical way in a fiat universe with cosmologi- 
cal constant. For t he details of the calculation, see 
lHallman et "all (j2007l ). 

The RMS amplitude of the density fluctuations as a 
function of mass and redshift, smoothed by a spherically 
symmetric window function with comoving radius R, can 
be computed from the matter power spectrum using: 



a 2 (M,z) 



^ P{k ,z)\Wn(k)\ 2 (6) 



where W#(/c) is the Fourier transform of the real-space 
top hat smoothing function: 



Wfl(fc) 



k 3 R 3 



[sin(kR) - kRcos(kR)} 



(7) 



We choose the value of gm for our lower limit from 
the lower mass limit at z=l, which is Muma = 1O 14 M . 
Simply put, we use a constant gm limit, and take all 
clusters in the simulation above that lower mass limit at 
each redshift. The result is shown in Figure [71 

The result clearly changes the sense of the redshift 
trend from the mass-limited sample. Now we see an in- 
crease in the fraction of clusters with cold fronts as a 
function of time. This selection should allow us to isolate 
evolutionary effects, as we are sampling the same part of 
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Fig. 6. — Mass limits derived from the gm~ and flux-limited 
cases described in the text. In each case, the sample of clusters 
used in the analysis are those above the specified mass limit at 
each redshift. 

the mass function at each redshift. The fraction gener- 
ally increases as we move forward in time, which is not 
particularly surprising, since the result of our sample se- 
lection is that the lower mass limit goes up with time. As 
we have shown, the fraction of clusters with cold fronts 
is a strong function of mass. What is puzzling is the dip 
in the fraction at z=0.2. It is not entirely clear how to 
explain the appearance of this feature, though it is possi- 
ble that a varying lower mass limit may create arbitrary 
features. 

The second additional selection mimics an observa- 
tional selection by a flux limit. Since this simulation 
has only adiabatic physics, and as such does not predict 
the X-ray luminosity accurately, we do not use the total 
X-ray emission calculated from the physical properties 
of the simulated clusters. Instead, we take a simple ap- 
proach, using observationally determined X-ray scaling 
relations to create a mass limit as a function of redshift 
for our sample. 

We take the X-ray scaling relation fit by iMantz et al.l 
(2009), a simple power law 

l(m) = p l m + (3[ m m, (8) 

where 

1 = 10910 Uz)iX g5 -0 ' (9) 

and 

m = l09w {-Wn^)- (10) 

E(z) is defined in the usual way in a flat ACDM universe, 

E{z) = {£l m {l + zf + £l A ) 1 / 2 , (11) 

and /ffi 71 and are the fitted parameters. IMantz et al.l 
(|2009l ) find those parameters to be P l m = 0.82±0.11, and 
pirn _ i 29±0.07 by fitting to the X-ray luminosity func- 
tion. We use the flux limit from the ROSAT-ESO Flux - 
Limited X-ray sample f REFLEX: [Bohringer et al.ll200lf ). 
In this case, we are not interested in any particular flux 
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Fig. 7. — Trend in redshift of clusters with the indicated L c f for 
each line. Dotted line indicates L c f > 78/i -1 fcpc, solid line is for 
our fiducial limiting L c f > 156h~ 1 kpc, dashed line is for L c f > 
2>12h~ 1 kpc. This plot uses a lower mass limit in each redshift bin 
drawn from a constant gm limit, highlighting evolutionary trends 
in the same set of clusters through time. Included are the trends 
with the same limiting length of cold front pixels as in previous 
Figure 2] Gray regions indicate la Poisson error bars. 

limit, we simply want to show the effect of sample se- 
lection on the cold front statistics. The REFLEX flux 
limit is 2.0 x 10~ 12 er gs~ 1 cm~ 2 . At each redshift from 
the simulation, we calculate the X-ray luminosity appro- 
priate to the flux limit, then calculate the cluster mass 
assoc iated to that lumin osity using the scaling relation 
from iMantz et al.l (|2QQ9h . For the z=0 clusters, we use 
£=0.05 to calculate the value of Dl and the flux. We 
compare the lower mass limits as a function of redshift 
in the gm~ and flux-limited cases in Figure El 

The resulting statistics of cold front incidence are 
shown in Figure [51 Note that this plot only shows frac- 
tions of clusters with cold fronts from z=0.5 to z=0, as 
at this flux limit, we have no clusters in our simulation 
above the appropriate mass at higher redshift s. In any 
case, it is clear that a flux limited selection results in 
a very different picture of the incidence of cold fronts 
when compared either to the mass-limited or dM-limited 
cases. The flux limit, of course, preferentially selects 
higher mass objects at higher z. As we have seen in 
our analysis of the incidence of cold fronts as a function 
of mass, it is a strong function of mass, so this result 
is not surprising. In a flux limited sample, we should 
expect the highest fraction of cold front clusters at high 
redshifts (since we select the highest mass clusters there), 
and it should decrease monotonically as we move to lower 
redshift. In other words, given that the flux limit prefers 
high mass clusters at high redshift, and we know that 
there are more cold fronts in high mass clusters at all 
redshifts from our analysis, we expect the highest frac- 
tion of cold front clusters at high z in a flux-limited sam- 
ple. In this selection, the total fraction of clusters with 
L c f > 156/i — 1 kpc in the full sample is ~15%, lower than 
typically found in observed samples. 

4.3. Properties of the Cold Fronts 
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Fig. 8. — Trend in redshift of clusters with the indicated L c f for 
each line. This plot uses a lower mass limit in each redshift bin 
drawn from a flux limit (from the REFLE X clus ter sample) and the 
X-ray scaling relations from Mantz et al. ( 2009). Each cluster mass 
is converted to a luminosity from the scaling relations, and then a 
flux given its luminosity distance. Included are the trends with the 
same limiting length of cold front pixels as in previous Figures 3] 
and [3 This includes our fiducial limiting L c f > 156h~ 1 kpc (solid 
green line), as well as limits of half the fiducial length (dotted, 
blue line) and twice the fiducial length (dashed, red line). Yellow, 
gray, and orange regions are la Poisson error bars for the L c f > 
78h~ 1 kpc line, L c f > 156h~ 1 kpc line and the L c f > 312h~ 1 kpc 
line respectively. 

Here, we study the properties of the identified simu- 
lated cold fronts - in particular the jump in temperature 
and surface brightness across the features. We expect 
from previous work looking at the sho ck population (e.g., 
iRvu et al.ll2QQl ISkillman et al.ll2008D . that the distribu- 
tion of number of pixels as a function of ratio of T s i and 
0.3-8 keV X-ray surface brightness, should be roughly a 
power law. As with shocks, in the formation of large-scale 
structure, we expect more small-ratio jump features than 
large-ratio ones, meaning that low velocity and temper- 
ature contrasts tend to be more plentiful than high con- 
trasts. With our simple cold front identifier, we indeed 
find that the temperature jump frequency for all the clus- 
ters from z=l to z=0 follow a rough power law with a 
break, as shown in Figure The jumps in surface bright- 
ness follow a less simple distribution, flat at low jump 
ratios and rolling off at large values as shown in Figure 
[TOl A population of pressure continuous features in the 
simulations should produce a distribution of temperature 
and density jumps with similar shapes. Though here we 
show the temperature and surface brightness jumps, un- 
der the assumption that the soft X-ray emissivity is only 
weakly dependent on temperature, the surface brightness 
and density jumps have very similar shapes (though den- 
sity is proportional to roughly the square root of X-ray 
emissivity). It is unclear why the shape of the tempera- 
ture jump distribution does not qualitatively match the 
shape of the surface brightness jumps. 

Additionally, we can check the number of identified 
cold front pixels per simulated cluster image. We have 
already described this measurement as equivalent to the 
comoving total length (L c f) of cold front features in the 
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Fig. 9. — Frequency of cold fronts as a function of ratio of T s i 
ratio from one side to the other identified in all cluster images 
(inside rsoo) from 1.0> z >0.0. 
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Fig. 10. — Frequency of cold fronts as a function of ratio of 0.3- 
8.0 keV X-ray surface brightness from one side to the other in all 
simulated cluster projections (inside rsoo in the range 1.0> z >0.0. 

full image, though we make no effort here to determine 
if the features are contiguous. For now, this calculation 
gives us a sense of the total extent of cold fronts in all 
the simulated cluster images. Figure [TT1 shows the result 
of this calculation, a probability distribution function for 
the total comoving length of cold fronts in the synthetic 
cluster maps at three redshift s. Interestingly, the distri- 
butions at three redshifts (1.0, 0.5 and 0.0) overlap quite 
well. The normalization is adjusted by the total num- 
ber of objects with cold fronts identified, but the relative 
number of objects with a specific measured length of cold 
fronts is more or less constant across redshift. The high 
L c f end of the distribution grows longer at later times, 
presumably because the objects progressively get larger, 
and therefore have larger values of rsoo, thus more area 
within which to find cold fronts. 
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Fig. 11. — Normalized probability distribution function for the 
comoving length of cold fronts (L c f) in all cluster projections at 3 
redshifts as indicated in the legend. The distribution of lengths is 
fairly constant with redshift, though the high L c f tail lengthens at 
lower redshift, due to the increased size of the largest clusters. 

Since cold front jumps are observed to be continuous in 
pressure, we check, in projection, whether these features 
are consistent with observations. If we take the projected 
surface brightness jump and use it to estimate (roughly) 
the projected density jump, and compare that to the cor- 
responding projected temperature jump, we can check 
to see if pressure continuity is met. The X-ray surface 
brightness in the 0.3-8.0 keV band is only very weakly 
dependent on temperature, and is effectively representa- 
tive of the projected value of n 2 e . For one of our simulated 
cold front clusters, we show the properties of a cold front 
in Figure [T21 We use the square root of the X-ray sur- 
face brightness jump as a proxy for the projected density, 
and multiply it by the value of T s i to form a projected 
pressure map. Figure [12] shows the obvious difference 
between the hot feature (a shock) and the cold feature 
behind it (a cold front). In the shock case, the jump in 
temperature corresponds to a pressure jump, and in the 
cold front case, the temperature jump is effectively con- 
tinuous in pressure. Note also the opposite direction of 
the jump in surface brightness and temperature in the 
cold front. These jumps are typical of the shocks and 
cold fronts in the simulated clusters. 

4.4. Cold Fronts in Individual Clusters 

The cold fronts identified in the images of the simulated 
clusters have properties consistent with those observed 
in real galaxy clusters. The jumps in temperature are 
typically factors of 1.5-3, with a corresponding inverse 
change in the X-ray surface brightness, such that the 
features are continuous in pressure. Cold fronts in the 
simulated images are typically associated with shocks, 
as is clear in Figure [12] While we have not specifically 
looked at the observability of these features in the X-ray, 
it is clear from our data that the surface brightness in 
the hot shocked region is significantly lower than that in 
the cluster core. It is possible that many of these shocks 
would not be detected in X-rays for a typical exposure 
time, though the cold front would be more easily de- 



tectable. In future work, we will explore the statistics of 
cold fronts given the X-ray observability using more so- 
phisticated synthetic observations including instrument 
response and X-ray backgrounds. 

It is also worth exploring the recent history of the sim- 
ulated clusters with the most obvious cold front features. 
Simulations show that cold fronts shou ld be a good diag- 
nostic for recent merger activ ity (e.g. jMathis et all 2005; 
IZuHone fc Markevitchll2009| ) , both in the case of sloshing 
type cold fronts and the more extended merger type cold 
fronts. A multipanel image of some of the most extended 
cold fronts in our simulated sample at z=0 are shown 
in Figure [13] The left column of this image shows the 
projected spectroscopic-like temperature (with identified 
cold fronts shown as white contours), and the right col- 
umn is the X-ray surface brightness. In these cases it is 
very clear that the cold front finder has done a sufficiently 
accurate job of characterizing the cold front regions. We 
note that these cold fronts are all of the merger type, and 
not the sloshing type. We do not see the sloshing type 
cold fronts in these adiabatic clusters, and do not expect 
them given the lack of radiative cooling and steep central 
entropy gradients. 

For three of these obvious cold front clusters, we have 
taken snapshots of their recent history to illustrate the 
process by which these extended cold fronts are formed in 
the cosmological simulations. Recall that for this large 
set of simulated clusters, the physics in the simulation 
is purely adiabatic, meaning that there is no radiative 
cooling. Therefore any regions where cold and hot gas are 
in close proximity result either from shocks or adiabatic 
effects. The histories of these three clusters are shown 
in Figures 1141151 and [16] In each case, we see that the 
formation of the cold fronts is a result of a recent merger. 

In Figures HH and [16j we show two classic roughly 
plane-of-the-sky mergers, where in the earlier timesteps, 
the halos of the merger are clearly approaching each 
other, as evidenced by the shock features between them. 
At the later epoch, we see the subclusters after core pas- 
sage, when the cold fronts have developed. We also show 
the dark matter contours, and there appears to be some 
mild separation of the dark matter centroids from the 
gas peaks, though we have not quantified this result in 
this work. In Figure [15j we see an array of cold fronts 
throughout the history of the cluster, but in this case the 
dynamics are more complicated. In the far left panels, 
we see a situation where there has been an initial core 
passage of two subclusters. The second panels show the 
approach of the two objects as they fall back together, 
with an obvious shock between, while at the same time 
a third subcluster approaches from the lower left. In the 
third panel, the initial merger has had a second core pas- 
sage, while the extra merging subcluster is now driving 
a shock as it wraps around the subcluster in the upper 
left. At the final stage, there are two very obvious cold 
fronts, but deducing the merging history from the final 
snapshot could be quite challenging. 

Cold fronts in the simulated clusters result from a va- 
riety of dynamical scenarios, ranging from simple binary 
mergers to multiple simultaneous mergers. The appear- 
ance of cold fronts and shocks is however a generic result 
of all kinds of mergers in these systems. 

4.5. Cold-Front Cluster Morphology 
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Fig. 12. — Left to right: Projected images of T s i, 0.3-8.0 keV X-ray surface brightness, and pressure for a cold front cluster with M = 
2x10 15 Mq at z=0 in the SFLC simulation. Also shown is the profile of T s i, X-ray surface brightness and projected pressure across the 
shock and cold front in the lower panels. Location of the profile pixels is shown in the T s i image, where pixel zero is at the bottom of the 
line, and pixel 36 is at the top in the image. Pixel scale is 15.6h _1 kpc. Note the very hot shocked gas and the strong drop in temperature 
at the same location where the X-ray surface brightness goes up sharply behind the shock. Also note that at the shock, there is a strong 
pressure jump, while at the cold front (where the temperature drops and the surface brightness peaks) the pressure is roughly constant. 



Recent studies of well-known cold front clusters iden- 
tify a possible relationship between the presence of X- 
ray cold fronts an d the merger state of the clusters 
([Owers et al.ll2009[ ). There is also some indication that 
in some fraction of clusters the presence of the cold 
front may be the on ly indicator of a recent merger 
(|Ghizzardi et al.l feOlO). To test these ideas, we have un- 
dertaken a study of the X-ray morphology of the clusters 
in our simulated samples via the measurement of power 
ratios and centroid shifts of the X-ray images. These 
quantitative measures are an indicator of whether the X- 
ray morphologies are obviously disturbed or regular, and 
are a proxy for de ciding whether a cluster is dynamically 
activ e or relaxed (|Jeltema et al J [20081 : IVentimiglia et aLl 
l2008h . 

In short, the power ratios are the multipole moments of 
the X-ray surface brightness map from a cir cular aperture 
centered on the cluster X-ray ce ntroid (see lBuote fc Tsai] 
1995. ll996tlJeltema et al.ll2QQ5h . As described elsewhere, 
this method is related to the multipole expansion of 
the two-dimensional gravitati onal potential . This set of 
equations is reproduced from lJeltema et all ([2008) . The 
multipole expansion of the two-dimensional gravitational 
potential is 
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where x' = (R', <//) and E is the surface mass density. 
In the case of X-ray studies, X-ray surface brightness 
replaces surface mass density in the calculation of the 
power ratios. 

The powers are formed by integrating the magnitude 
of \I/ m , the rath term in the multipole expansion of the 
potential given in equation (12), over a circle of radius 
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The higher multipole moments are sensitive to the 
amount of substructure in the cluster. In most cases, 
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Fig. 13. — Images of the simulated galaxy clusters with the 
most extended cold front structures. Left column is projected 
spectroscopic-like temperature with cold fronts identified using the 
automated method described in the text overlaid as white con- 
tours, second column is 0.3-8.0 keV X-ray surface brightness. Im- 
age scale is ~3 Mpc. Colorbar indicates temperature scale from 
left column in keV. Masses of the clusters from top to bottom are: 
1.4xl0 15 M©, 8.5xl0 14 M©, 4.8xl0 14 M©, 1.98xlO 15 M . 

as in this work, we normalize the multipole moments P2, 
P3 and P4 to Pq to factor out the overall brightness of 
the cluster. P2, the second multipole moment is sensi- 
tive to deviations from circularity, P3 to deviations from 
mirror symmetry, and P4 is similar to P2, but sensitive 
to smaller scale asymmetry. Each of the moments has a 
component of radial decline in the cluster profile which 
contributes to its value, such that substructure at large 
radius leads to larger power ratios. Perfectly round clus- 
ters (irrespective of the shape of the radial profile) would 
give all zero power ratios. While each of the power ratios 
is sensitive to slightly different asymmetries, it is impor- 
tant to note that they are strongly correlated with one 
another. That means that clusters that are disturbed via 
merging activity typically show deviations from circular- 
ity, as well as from mirror symmetry and small scale sub- 
structure simultaneously. In this work, we use the values 
of the power ratios inside a radius of rsoo, a radius which 
is typically observable in X-rays for a large set of galaxy 
clusters. We note that the value of the power ratio se- 
lected is somewhat dependent on the outer radius used, 



11 




0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 



Fig. 14. — Time sequence of development of a cold front in a 
1.4 x 1O 15 M cluster (same as the cluster in the top panels of Figure 
[13]). Left panels are at 2=0.1, right are at z=0. Upper panels show 
projected spectroscopic-like temperature in color and 0.3-8.0 keV 
X-ray surface brightness in contours. Color range in T s i goes from 
3 keV (blue) to 10 keV (white). Bottom panel is X-ray surface 
brightness in color and projected dark matter density contours. 
Colorbar on this image shows range of X-ray surface brightness 
in er gs / s / cm 2 /pixel at the face of the projected simulation box 
(not corrected for distance of observer) to allow a comparison be- 
tween the two images. Image scale is 2.0 Mpc. Note the obvious 
merger at z=0.1, which creates the cold fronts after the cores of 
the subclusters pass through each other. 



TABLE 1 

Median Power Ratios and Centroid Shifts 



Measure 


Cold Front Clusters 


Non-Cold Front Clusters 


P2/P0 


1.75e-6 


5.71e-7 


P3/P0 


2.25e-8 


4.02e-9 


P4/P0 


6.50e-9 


8.57e-10 


Centroid Shift 


0.0146 


0.0080 



and that merging clusters with disturbances on larger 
scale s than the outer radiu s may be missed when using 
r 500 (jTakizawa et al.ll2010D . 

A second type of quantitative mea sure we employ in 
this study is the centroid shift (see Mohr et al.l 119931 : 
lO'Hara et all 120061 : iPoole et all 120061 : iMaughanl f2QQ7h . 
Centroids of the X-ray surface brightness are measured 
successively in apertures of increasing size. The cen- 
troid shift is the variati on of the centroid w ith radius 
(see above references and lJeltema et al.ll2QQ8h . and is an 
indicator of deviations from dynamical equilibrium. Our 
main interest here is in the change of the mean centroid 
shifts from the non-cold front clusters to the cold front 
clusters. 

We find that clusters identified as having cold fronts 
have systematically higher values for the power ratios 
and centroid shifts than those that do not have cold 
fronts. This result is statistically significant at > 99% 
when a Kolmogorov-Smirnov test is applied. Table [T] 
shows the median values for the power ratios and cen- 
troid shifts for clusters identified as cold front clusters 
and those without cold fronts. Figures [T7] and [18] show 
the distribution of P3 / Pq and centroid shifts in the clus- 






Fig. 15. — Time sequence for a cold front developing in a 1.98 xlO 15 M0 cluster (same as the cluster in the bottom panels of Figure 
I13|) . Left to right panels are forward in time, from z=0.3, z=0.2, z=0.1 and z=0.0. Upper and lower panels show the same quantities as 
in Figure [14] Image scale is 4.0 Mpc, colorbar is T s i in Kelvin for the upper panels. Note the obvious cold fronts in the earliest timestep, 
followed by a recollapse and shock, then an additional subgroup merges at the same time in the third timestep, causing the development 
of multiple cold fronts simultaneously. 




Fig. 16. — Time sequence for a 4.8xlO 14 M cluster same as 
the cluster in the second from bottom panels of Figure [13]) , right 
panels are z=0, left are z=0.1. Image quantities are the same as 
Figures fl4] and [T5l Colorbar is T s i in Kelvin for the upper panels. 
Image scale is 2.0 Mpc. Note the similarities to Figure [14] with 
an obvious plane of the sky merger in the first timestep and the 
classic post core-crossing cold fronts in the second timestep. 



ter sample for both cold front and non-cold front clusters 
in the sample at z=0. The other power ratios have very 
similar distributions. 

This result holds for all clusters with even a single pixel 
identified as a cold front, but it becomes monotonically 
more significant when we increase the threshold for num- 
ber of pixels identified as cold fronts. This measure is 
equivalent to total length of cold front features. Figure 
n~UI shows a histogram of the P3 / Pq power ratio for the 
clusters in samples with increasing numbers of cold front 
pixels. As we raise the threshold for number of pixels 
identified as cold fronts, the morphology becomes sys- 
tematically more disturbed, as is expected if cold fronts 



0.10 

0.05 
0.00 



I I No Cold Fronts 
CI2 Cold Front Clusters 



log(P3/P0) 



Fig. 17. — Probability density function of clusters as a function 
of P3/P0 power ratio for each cluster image in the z=0 sample. 
Green dashed line is cold front clusters, blue solid line is non-cold 
front clusters. Cold front clusters are those that have a minimum 
comoving length of cold front pixels L c f > 156/i -1 kpc comoving. 



are related to merging events. This trend is repeated 
in the other power ratios, as well as the centroid shift 
measure. In Figure [19j the number of cold front pixels 
identified in the z=0 simulated clusters increases as the 
distribution shifts to the more disturbed (right) side of 
the distribution. As noted in earlier sections, L c f can be 
thought of roughly as a total extent of cold front features, 
where in these images, the pixel scale is 15.6/i _1 kpc. 

5. EFFECTS OF BARYONIC PHYSICS 
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Fig. 18. — Probability density function of clusters as a func- 
tion of centroid shift for each cluster image in the z=0 sample. 
Green dashed line is cold front clusters, blue solid line is non-cold 
front clusters. Cold front clusters are those that have a minimum 
comoving length of cold front pixels L c f > 156/i _1 kpc comoving. 
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Fig. 19. — Cumulative fraction of cluster images as a func- 
tion of Ps/Pq power ratio for clusters with increasing thresh- 
old for number of cold front pixels identified. Blue line repre- 
sents all cluster images with no cold front pixels identified. Red 
line represents 5 identified pixels (L c f > 78/i _1 kpc), green line 
is for 10 pixels (L c f > 156h~ 1 kpc), cyan line is for 20 pixels 
(L c f > 312h~ 1 kpc), purple line is for images with 50 or more 
cold front pixels (L c f > 780/i — 1 kpc). Note the shift in the median 
of the distribution to more disturbed (higher) P3 / Po values as the 
threshold for minimum number of cold front pixels is raised. 

Two other simulations are used to study the effect of 
additional non-gravitational physics on the properties of 
cold fronts. While not a convergence study, these results 
are presented to give a basic indication of the impact 
of additional baryonic physics on the result. We know 
observationally that in galaxy clusters the gas physics is 
not purely adiabatic. It is well-known that in the core 
regions of clusters, radiative cooling becomes important, 
in many clusters the estimated cooling times are signif- 



FlG. 20. — Fraction of cluster projections hosting at least 10 
pixels (L c f > 156h~ 1 kpc) identified by the automated cold front 
finder. Solid line is for the CC-OTH run described in the text, 
with radiative cooling and star formation plus feedback. Dashed 
line represents the SFLC simulation, and is replotted data from 
Figure [4] Dotted line is for the CC-Adia run, identical to CC- 
OTH, but with only adiabatic gas physics. Gray regions are la 
Poisson error bars. 

icantly shorter than a Hubble time. In addition, there 
are both heat and metals injected from cluster galaxies 
due to star formation, and there is also strong evidence 
of AGN activity. In this case, we explore the result of 
including additional physical processes on the incidence 
of cold fronts in the simulated clusters. 

Our main result is shown in Figure [20] showing the 
fraction of clusters with at least 10 cold front pixels (in 
this case the images have the same pixel scale across all 
simulations) as a function of cluster mass. The plot axes 
match Figure 2J and the result from that figure are over- 
plotted with the same result from the CC-Adia and CC- 
OTH simulations. Several things are obvious in this com- 
parison, the first is that these simulations probe different 
cluster mass ranges, since the SFLC represents a much 
larger physical volume than CC-Adia and CC-OTH. The 
second is that including radiative cooling and feedback 
results in significantly more clusters hosting cold fronts. 
This is not unexpected, as the cooling of cluster cores, as 
well as subclusters, results in more cold edges surrounded 
by hot gas. Figure [21] illustrates this effect, comparing 
the identical cluster in both the CC-Adia and CC-OTH 
runs. 

Figure [22] shows the trending of fraction of cold front 
clusters with redshift in our simulation with additional 
baryonic physics (CC-OTH). While the fraction has some 
variation with redshift, within the statistical errors of the 
sample the trend is roughly flat. This compares well with 
the SFLC adiabatic run, though there is a mild increase 
in the fraction of cold front clusters from high to low 
redshift in that simulation. In CC-OTH, the fraction 
goes from roughly 0.64 at z=0.9 to 0.54 at z=0. However, 
as noted, the Poisson errors are large (~ ±0.1 at all z) 
due to the smaller cluster sample in this simulation. 

Additionally we can compare the fraction of clusters 
with cold fronts from SFLC to CC-Adia, since the physics 
is identical, though the mass resolution is different. CC- 
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Fig. 21. — Comparison of T s i maps for the same cluster in two different simulations, the left panel is from the simulation designated 
CC-OTH and the right is from CC-Adia. T s i is the image quantity, the contours are synthetic 0.3-8.0 keV X-ray surface brightness contours 
with the same limits in each map. The 6.1xlO 14 M0 cluster from the CC-OTH (simulation with radiative cooling, star formation and 
thermal/metal feedback) in the left panel, and the same cluster in the adiabatic simulation is shown in the right panel. Image scale is 
roughly 3/i -1 Mpc square for each panel. Note the significant difference in the temperature maps, with the cold subclumps very obvious in 
the CC-OTH cluster map. Note that the cold subclumps create peaks in the X-ray contours. 

We have shown the incidence of X-ray cold fronts in 
a large sample of simulated galaxy clusters. In this ini- 
tial study, to maximize the number of clusters in the 
sample at useful grid resolution, we gather the statistical 
results from an adiabatic simulation. We identify cold 
fronts using projected maps of spectroscopic-like tem- 
perature (T s i ), and X-ray surface brightness generated 
with Cloudy. We show the effect of additional baryonic 
physics using two smaller volume simulations, identical 
in initial conditions, one using adiabatic physics, and the 
other including the effects of radiative cooling, star for- 
mation and metal/thermal feedback. The main results 
are summarized here. 

• The trends of simulated galaxy clusters with cold 
front features with redshift depends significantly 
on cluster selection. When using a mass-limited 
sample the trend out to z=1.0 is weak, irrespective 
of the threshold of number of identified cold front 
pixels. This weak trending occurs only when the 
sample is mass-limited, however. When limiting 
the sample by the rarity of peaks using a constant 
gm , the incidence of cold fronts shows a decreasing 
trend with redshift, and when using a flux limit 
derived from X-ray scaling relations, the trend is 
increasing with redshift. 

• There is in the simulated cold front clusters a 
strong trend with cluster mass, the more massive 
clusters being much more likely to host a cold front. 
In the adiabatic sample, the fraction of clusters 
hosting cold fronts with L c f > 156h~ 1 kpc with 
M> 7.5 x 1O 14 M is 40-50% depending on the red- 
shift of the sample. In the simulation with radiative 
cooling, star formation and thermal feedback, that 




Fig. 22. Fraction of clusters with M> 10 14 in the CC-OTH 
simulation with cold front length L c f > 156h~ 1 kpc as a function 
of redshift. Gray region is la Poisson error bars. 

Adia has a lower incidence of cold front clusters than does 
SFLC, which is somewhat surprising. However, given 
the number of changes between SFLC and CC-Adia, it 
is premature to make strong conclusions. For example 
the much reduced box size in our smaller simulations 
(from 512/i _1 Mpc to 128/i _1 Mpc cubic volume) leads to 
a radical reduction in large-scale power in the simulation 
box. What has not changed between any of these runs, 
however, is the increasing trend of fraction of cold front 
clusters with mass. 

6. DISCUSSION AND SUMMARY 
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fraction can be higher than 80% at the high mass 
end. 

• The fraction of clusters hosting cold fronts in the 
adiabatic simulations is lower than in observed 
samples, which give fractions anywhere from 40- 
80% depending on the sample and the redshift cov- 
erage. We know given our method that we proba- 
bly underestimate the total number of cold fronts, 
so this is not surprising. In our flux-limited sam- 
ple from the simulations, we find 15% of cluster 
with L c f > 156/i _1 kpc. We show that including 
radiative cooling and star formation in the simula- 
tions will increase this fraction significantly with- 
out changing the general trends with mass and 
redshift. Additionally as noted above, the high- 
est mass clusters in our sample show a higher cold 
front fraction, around 40-50%. 

• Cold fronts in simulated clusters are almost uni- 
versally the result of mergers of various sizes, and 
are almost always associated with shocks. This re- 
sult is still qualitative, but we hope to make more 
quantitative claims in subsequent work. 

• The association of cold front clusters with mergers 
is confirmed in a statistical sense using the quan- 
titative morphological measures for the X-ray sur- 
face brightness maps. Clusters with cold fronts are 
preferentially more disturbed when measured using 
power ratios and centroid shifts. The significance 
of this result increases if we select clusters with 
higher numbers of identified cold front pixels. 

• There is a clear impact on the incidence of cold 
fronts due to changes in baryonic physics of the 
simulation. When including radiative cooling and 
star formation and feedback to the simulations, the 
incidence of cold fronts goes up strongly to greater 
than 50% of clusters at all masses. While the nor- 
malization of the number of cold fronts in simu- 
lated clusters is strongly affected by the baryonic 
physics, the trending with mass and redshift is not. 

The most obvious improvements that can be made to 
this analysis are a full resolution and convergence study 
with variations in baryonic physics, and an analysis of 
the statistics of cold fronts expected with more realis- 
tic X-ray images and temperature maps including X-ray 
instrumental effects and backgrounds. 
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